Известия Томского политехнического университета. 2006. Т. 309. № 1 


УДК 539.125.5 

КРИТИЧЕСКИЕ НЕЙТРОННО-ФИЗИЧЕСКИЕ ПАРАМЕТРЫ УРАН-ТОРИЕВЫХ 
И ПЛУТОНИЙ-ТОРИЕВЫХ СПЛАВОВ 

В. И. Бойко, П.М. Гаврилов*, И. В. Шаманин, М.Г. Герасим*, В.Н. Нестеров 

Томский политехнический университет, г. Томск 
**ФГУП Сибирский химический комбинат, г. Северск 
Е-таіІ: $Ііатапіп@к21.р)іШри.ги 

Анализируется возможность безопасного хранения сплавов сырьевого нуклида V т 232 с основными нечетно-четными нуклидами 
ІІ Ш и Ри 239 . Получены соотношения для определения предельно допустимых значений концентрации ядер урана и плутония в 
сплавах, приведены результаты нейтронно-физических расчетов. 


Введение 

В результате реализации уран-ториевого топлив- 
ного цикла существенно улучшаются технико-эко- 
номические показатели, повышается безопасность 
АЭС с реакторами на тепловых и быстрых нейтро- 
нах, закладываются основы ядерной энергетики бу- 
дущего. Уран-ториевый цикл позволяет снизить на- 
копление плутония, трансурановых элементов и ак- 
тиноидов, играющих негативную роль при решении 
вопросов захоронения радиоактивных отходов [1]. 

Концепция открытого уран-торий-плутониево- 
го топливного цикла и реакторных установок на 
его основе не противоречит имеющимся достиже- 
ниям и опыту развития ядерной энергетики. На- 
против, предлагаемая концепция расширяет перс- 
пективы развития уже освоенных отечественной 
промышленностью традиционных типов реакто- 
ров (прежде всего легководных реакторов и реакто- 
ров на быстрых нейтронах) [2-4]. 

В связи с перспективой использования оружей- 
ных урана и плутония в ядерно-энергетических 
установках актуальным является вопрос о ядерной 
безопасности хранения и изготовления биметалли- 
ческих соединений: (яІР 35 тТй 232 ) и (яРи 239 тТЪ 232 ). 
Возникает вопрос - какова доля урана или плуто- 
ния в смеси с торием, которая допустима правила- 
ми ядерной безопасности? Будет ли данная доля 
приемлема с точки зрения технического обоснова- 
ния эффективности технологии изготовления би- 
нарных сплавов? Поиск ответов на данные вопро- 
сы является задачей данной работы. 

Необходимое условие самопроизвольного деления 

смеси основных делящихся нечетно-четных нукли- 
дов с сырьевым четно-четным нуклидом ТН 232 

Под основными нечетно-четными нуклидами в 
работе подразумеваются ІР 35 и Ри 239 , являющиеся 
базовыми компонентами ядерного топлива. Суще- 
ство вопроса составляют рассмотрение замедления 
нейтронов, образующихся при делении ядер урана 
и плутония, и определение условий, необходимых 
для возникновения цепной реакции деления урана 
и плутония, смешанного с торием. 

Для возможности развития цепной реакции де- 
ления основных нуклидов необходимо, чтобы об- 


разующиеся при делении ядер урана и плутония 
быстрые нейтроны успевали с достаточной вероят- 
ностью вызвать следующий акт деления до того, 
как они покинут объем биметаллической смеси. 
Если же образующиеся при делении нейтроны за 
“время жизни” успевают замедлиться до энергии, 
для которой сечение деления максимально (тепло- 
вые нейтроны), то возможность возникновения 
неуправляемой цепной реакции существенно воз- 
растет. Ситуация упрощается тем, что нейтроны за- 
медляются в среде с большой атомной массой ядер, 
но осложняется тем, что “сырьевой” нуклид ТЬ 232 
эффективно делится нейтронами высоких энергий. 

Ответ на вопрос о “цепном распаде основного 
изотопа урана” найден и представлен в классиче- 
ской работе [5]. Методология данной работы ис- 
пользована при получении результатов, о которых 
ниже пойдет речь. Отличия состоят в том, что в на- 
стоящей работе рассматриваются смеси основных 
нуклидов с нуклидом, делящимся быстрыми ней- 
тронами. Деление быстрыми нейтронами учитыва- 
ется для всех ядер, присутствующих в биметалли- 
ческой смеси. 

Итак, необходимо найти вероятность со того, 
что образующиеся с начальной энергией Е п нейтро- 
ны без поглощения, в результате которого произой- 
дет деление, будут замедлены до энергии Е ь меньше 
которой нейтроны не могут иметь энергию, так как 
вызовут деление либо будут поглощены без деле- 
ния. В "тяжелой" среде значение Е к несоизмеримо 
выше привычной тепловой энергии. Это связано с 
тем, что, хотя микроскопическое сечение деления 
быстрыми нейтронами основных нуклидов очень 
мало по сравнению с таковыми для тепловых ней- 
тронов, из-за высокой (по сравнению с реакторным 
топливом) концентрации делящихся ядер макро- 
скопическое сечение оказывается значительным. В 
результате, быстрый нейтрон “не успевает” замед- 
лится до тепловой энергии, но в процессе много- 
кратных взаимодействий с малыми передачами 
энергии все же вызывает деление. В условиях пре- 
дельно большой массы смеси основных нечетно- 
четных нуклидов с сырьевым четно-четным нукли- 
дом Тй 232 необходимое условие инициации цепной 
реакции определяется неравенством: 

ѵ(1 — со) > 1, (1) 
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где ѵ - число нейтронов, образующихся на один 
акт захвата “первичного” нейтрона с энергией в 
области Е 0 ...Е к . 

Будем “отслеживать” историю нейтрона с уче- 
том упругого и неупругого рассеяния на ядрах до 
того момента, пока он не вызовет деление или не 
будет поглощен без деления. Выполнение таких 
условий, как и в работе [5], является необходимым, 
но может быть недостаточным для возбуждения 
цепной самоподдерживающейся реакции деления. 

Для определения значения со необходимо совме- 
стно решить уравнение, описывающее диссипацию 
энергии нейтронов при рассеянии, и уравнение, опи- 
сывающее скорость уменьшения количества нейтро- 
нов из-за поглощения, в том числе с делением: 

&—Е 

" = < 2 > 

I 


где ѵ - скорость нейтронов; Е - их энергия; N - 
плотность; Е е/ . - макроскопическое сечение упруго- 
го рассеяния на ядре сорта /; 8 е1і - доля энергии, те- 
ряемой при этом; Л іп ., 5 ІЩ - то же для неупругого рас- 
сеяния; Ё а , - макроскопическое сечение поглоще- 
ния /-ого сорта ядра (суммирование проводится по 
всем видам присутствующих в смеси ядер). Допу- 
ская, что о е1 , ст„„ сх„, 5 е1 , 5 ІП не зависят от энергии ней- 
тронов в рассматриваемом энергетическом интер- 
вале, систему ур. (2) можно преобразовать к виду: 


СІЕ П^'2 Г' 

-Ж~ Е с - 


^- = ЫЕ Ѵ2 С 2 , 


( 3 ) 


где функции С, и С 2 определяются только составом 
смеси нуклидов и могут корректироваться при из- 
менении значений Е п и Е к в соответствии с измене- 
нием значений сечения взаимодействия. С, и С 2 
определяются соотношениями: 

С 1 = -727^Х( Е 4^+2,,5 1 , 1 ), 


С 2 =-Ѵ27^Х, (4) 


где т - масса нейтрона. 

Решение системы (3) имеет вид: 



(5) 


где ІѴ/ЛГо =со. 

Величина ѵ в условии (1) определяется соотно- 
шением: 


V 


I У/° л рі 
IX А 


(6) 


где р; - концентрация /-го сорта ядра; а г и о а - ми- 
кроскопические сечения деления и поглощения, со- 
ответственно; Ѵ( - число нейтронов на один акт де- 


ления. Ядерно-физические характеристики необхо- 
димо задавать в энергетическом интервале Е^.-.Е^. 

Если допустить, что упругое рассеяние изотроп- 
но, а средний косинус угла рассеяния нейтрона в си- 
стеме центра инерции равен 0, то в соответствии с [6]: 


5„ = 


2 А 


(А + іу 


<5,„ = ■ 


2 А 


( А + 1 ) 2 
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> (7) 


где А - атомная масса ядра мишени; 0 - энергия 
возбуждения составного ядра при неупругом рассе- 
янии. Вычисления по соотношениям (1, 4-7) по- 
зволяют определить способность той или иной 
мультиплицирующей системы (среды) к возбужде- 
нию в ней цепной реакции деления. 


Определение "времени жизни" и конечной энергии 
нейтрона в тяжелой размножающей среде 

Для того, чтобы определить способность той 
или иной мультиплицирующей системы (среды) к 
возбуждению в ней цепной реакции деления, 
необходимо найти вероятность со , которая опреде- 
ляется соотношением: 



Итак, необходимо определить энергию Е к , до 
которой успевает замедлиться нейтрон от некото- 
рой начальной энергии Е 0 в результате упругих и 
неупругих столкновений с ядрами среды, состоя- 
щей из смеси изотопов Тіг 232 - И 235 - Ри 239 , а также 
“время жизни” г нейтрона в этой среде. Значение 
Е к определим используя численный эксперимент 
(метод Монте-Карло). 

Связь между энергией нейтрона до рассеяния Е' 
и после рассеяния Е имеет вид: 


Е = Е 


1 -- 


2^4 


( А + 1 У 


1 -А+10-+ ц 

2 Е' Ис 


А + 1 0 
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А Е' 


где /л с - косинус угла рассеяния в системе центра 
инерции. 

Возбуждение составного ядра в результате не- 
упругого столкновения с нейтроном возможно 
лишь при 


Е'> ^ 


(А + 1) 
А 


с вероятностью Х, л /Х 5 . 

Рассеяние нейтрона на тяжелых ядрах изотроп- 
но при энергии нейтрона меньше 100 кэВ. При 
больших энергиях угловая диаграмма направлен- 
ности рассеяния нейтронов вытягивается вперед, 
что означает, что доля энергии, передаваемая ядру 
нейтроном при упругом рассеянии, уменьшается. 
Увеличение значения доли энергии нейтрона, пе- 
редаваемой ядру отдачи при рассеянии, эквива- 
лентно уменьшению его “времени жизни” в тяже- 
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лой среде, то есть ужесточению требования к 
мультиплицирующей системе в части возможности 
возбуждения цепной реакции деления. Воспользо- 
вавшись этим можно допустить, что упругое рассе- 
яние изотропно. 

Алгоритм метода статистических испытаний 
для определения “времени жизни” нейтрона в тя- 
желой среде отличается от стандартного [7] и име- 
ет следующую структуру: 

1) определение длины свободного пробега ней- 
трона; 

2) определение времени пролета нейтрона до 
столкновения с каким-либо ядром; 

3) выбор случайного числа и определение типа 
ядра, с которым взаимодействует нейтрон; 

4) выбор случайного числа и проверка условия по- 
глощения. Если оно выполняется, то расчет 
траектории заканчивается. Энергия, которой 
обладал нейтрон после последнего рассеяния, 
считается конечной, а сумма времен между все- 
ми взаимодействиями до момента поглощения 
считается “временем жизни” нейтрона в дан- 
ной среде; 

5) если условие пп. 4 не выполнено, считается, что 
нейтрон испытал рассеяние; определяется вид 
рассеяния (упругое или неупругое), а также кос- 
инус угла рассеяния; 

6) определение энергии нейтрона после рассеяния; 

7) повторение процедуры расчета, начиная с пп. 1. 

Проследив таким образом достаточное количе- 
ство историй нейтронов, можно вычислить среднее 
значение конечной энергии и “времени жизни” 
нейтронов для конкретной смеси изотопов. 

В качестве начальной энергии нейтрона в рас- 
четах была выбрана средняя энергия нейтронов 
спектра деления 2 МэВ, зависимость нейтронных 
сечений от энергии аппроксимирована 26-группо- 
вым приближением [8]. В качестве значений 0 вы- 
браны энергии возбуждения первого уровня. В 
частности, для ЕР 35 0= 14 кэВ, для Ри 239 0=7,85 кэВ, 
для ТЪ 232 0=50 кэВ [9]. Рассмотрены случаи попа- 
дания или образования в бесконечной гомогенной 
среде нейтронов с энергией 10 МэВ, значительно 
превышающей среднюю энергию нейтронов спек- 
тра деления. 

В теории расчета критических систем до сих 
пор существует ряд пробелов. Это, прежде всего, 
проблема оценки систематических и статистиче- 
ских погрешностей реакторных функционалов в 
размножающих средах без источников [10]. Метод 
Монте-Карло позволяет определить значение эф- 
фективного коэффициента размножения для фик- 
сированных размеров, геометрии и материального 
состава, а не наоборот. Выбор оптимального соче- 
тания последних без привлечения аналитических 
методов требует неоправданно большого количе- 
ства численных экспериментов. В этой связи была 
решена следующая задача. 


Критические размеры в средах: тІІ 235 (1-/п)ТИ 232 

и тРи 239 (1-т) ТН 232 

При расчетной оценке возможности проведе- 
ния работ с размножающими системами (РС) ис- 
пользуют количественный критерий ядерной безо- 
пасности [11]: 

К<К Ф <КІ = 

= +Х АД* +«[2>' 2 (ЛА, )] 1/2 } ^ К (9) 

і і 

где к ь к 2 - границы промежутка, внутри которого 
оценку условий по значению эффективного коэф- 
фициента размножения РС делают расчетным спо- 
собом; к*ф - предписываемая правилами ядерной бе- 
зопасности смещенная оценка к эф РС; к" 4 — рассчи- 
танное по специализированной вычислительной 
программе (в нашем случае - МС1ЧР, ЬА- 12625-М, 
Ѵешоп 4В) значение кф ХА ік зф - алгебраическая 
сумма поправок Д,^ к расчетному значению к р ф 
компенсирующих систематические погрешности 
расчета; а - квантиль, который принимает по пра- 
вилам ядерной безопасности значения около 3...4, 
характеризует надежность выполнения условия и 
определяет ширину интервала, накрывающего с за- 
данной вероятностью погрешность расчета; 
[Е5 2 (Л,Е^)] 1/2 - оценка среднеквадратичного откло- 
нения поправленного расчетного значения 
к Э ф=к р ,ф+Т,А,к Э ф, связанная со средними квадратиче- 
скими отклонениями <5(Л,А/,) каждой поправки А ,к эф . 
Сходные с приведенными в выражении (9) крите- 
рии ядерной безопасности сформулированы для ко- 
эффициентов по критической массе и критическим 
размерам. 

Случаю однородной среды с изотропным рассе- 
янием отвечает уравнение: 

0Ѵ(р(г,0) + Е ( (р(г,0) = 

= 0^+2 (10) 
Ап 

которое имеет решение на множестве К р в , 0=Ѵ/. Ѵ п . 
Принадлежность функции ср(г,0.) к К в означает, 
что она удовлетворяет граничному условию: 

<р(г 3 ,П) = 0, г е Г, ОеМ_(^). 

Задача об определении критических размеров 
состоит в том, что из множества сочетаний харак- 
терных размеров (ширинахдлинахвысота) выбира- 
ются элементы с заданными значениями мате- 
риальных параметров [12]. При этом в большин- 
стве интересных случаев получается либо пустое, 
либо одноэлементное множество. 

Пусть геометрия объема V фиксирована. Рас- 
смотрим вспомогательную задачу: 

0У^(г,0) + ^(г,0)=|^|^(г,0)й10, 


которая относится к разряду условно-критических. 
Эта задача всегда имеет положительное решение 
(к эф >0,\і/), которое может меняться лишь при изме- 
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нении геометрических характеристик V. Поэтому 
величину В г =\/к эф в (1 1) можно называть геометри- 
ческим параметром среды. 

Введем еще одну вспомогательную задачу: 


&Ѵ ф(г ,Сі) + ф(г ,СТ) = ^- 1 ф(г,С1)с1€1, 

4л 

Ѵу-Х г + X 

геѴ, В м = (12) 


где величину В и естественно интерпретировать как 
материальный параметр, так как она полностью 
определяется материальным составом среды и не 
зависит от ее геометрических характеристик. Ура- 
внение (12) получается из (10) путем введения но- 
вой переменной: г— >1., г. При этом новая про- 
странственная переменная измеряется в единицах 
длины полного свободного пробега нейтронов. 

При произвольном значении В м уравнение (12) 
не обладает положительным решением, но если 

В м =Ук эф , (13) 


где к эф — собственное число задачи (11), соответ- 
ствующее положительному решению, тоД 12) также 
имеет положительное решение [В м ,ф(г ,П)]. 

Критические размеры объема V при заданном 
значении с можно искать, решая задачу (11) для 
различных значений линейных размеров и выделяя 
затем среди них те, которые обеспечивают выпол- 
нение критического условия (13). 

Статистические испытания показывают, что 
энергия нейтронов лежит в одном энергетическом 
интервале 26-группового приближения. Это зна- 
чит, что, используя данное приближение можно 
считать сечения нейтронных реакций не завися- 
щими от энергии и постоянными в рассматривае- 
мом интервале энергии. 

Благодаря тому, что нейтронные сечения удо- 
влетворяют условиям приближения постоянных 
сечений, спектральная задача сводится интегриро- 
ванием по энергии к односкоростной задаче. Так 
как конечная энергия нейтрона лежит выше поро- 
говой энергии деления тория, то 


я( ѵа г + о-, ) + п тъ (ѵ ть о - у ТІ1 + с/ 11 ) 
м п(а г + <у, + ст с ) + я ть (р? + ст/ ь + ст™ ) : 
2 ,г = и(°7 +<т +<т)+ и^сгТЧст/Чсг™), 


где п, оу, ст„ ст„ ѵ - концентрация, микроскопиче- 
ские сечения деления, рассеяния (упругое + не- 
упругое), радиационного захвата, число вторичных 
нейтронов на один акт деления основного деляще- 
гося изотопа (ІР 35 либо Ри 239 ), соответственно; п ть , 
сг/\ сг/\ стД % - аналогично для ТВ 232 . 


Результаты расчетов 

На рис. 1-2 приведены расчетные зависимости 
“времени жизни” нейтронов и их конечной энер- 
гии в уран-ториевых и плутоний-ториевых сплавах 
от начальной энергии нейтронов и от содержания 
нечетно-четного нуклида. 


т, 10 ' с 



Рис. 1. Зависимости ". времени жизни " нейтрона в средах ви- 
да: 1) тР^пШ 32 и 2) тІ/^пТН 232 от концентрации 
плутония и урана при начальной энергии нейтрона: 
а)Ю МэВ ; 6) 2 МэВ ( точками обозначены результаты, 
полученные в численном эксперименте) 

“Время жизни” нейтронов с начальной энерги- 
ей 10 МэВ практически в 5 раз меньше, чем при на- 
чальной энергии 2 МэВ. Это обусловлено: 

- большей скоростью нейтронов в интервалах 
между актами рассеяния и, следовательно, 
большей частотой столкновений; 

- более высокими значениями микроскопиче- 
ских сечений в области высоких энергий. 

Зависимости, полученные в статистических 
экспериментах, аппроксимированы полиномами 
второго порядка. 

Для уран-ториевого сплава: 

1. При ^(,= 10 МэВ: 

г=4-10- 8 -2,5-10- 9 с+9-10- 11 с 2 ;^=8+0,08с-0,002с 2 . 

2. При Е 0 =2 МэВ: 

г=2 • 1 0 7 — 1 , 7 ■ 1 0 8 с+ 6 , 5 • 1 0 ^с 2 ; Е к =1+0, 04с-0, 001с 2 . 
Для плутоний-ториевого сплава: 

1. При ^0=10 МэВ: 

г=4- 10 8 -3- 10 9 с+ 1,5- 10-Ѵ; Я*=8+0, 1с-0, 004с 2 . 

2. При Е 0 =2 МэВ: 

г=2 ■ 1 0~ 7 -2 ■ 10 8 с+7- 1 0 _1 Ѵ; Е к = 1+0, 05 с-0, 001с 2 . 
Здесь с - концентрация нечетно-четного нукли- 
да в десятичных долях, г - в секундах, Е к - в МэВ. 

При подстановке полученных выражений для 
Е к в соотношение (8) определяется значение пара- 
метра (критерия), определяющего возможность 
инициации цепной реакции деления (1):Р=ѵ(1-йі). 
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Рис. 2. Зависимости конечной энергии нейтрона в средах 
вида: 1) тРі/^пТН 232 и 2) т^пТН 232 от концентрации 
плутония и урана при начальной энергии нейтрона: 
а) 10 МэВ ; 6) 2 МэВ 

На рис. 3 приведена зависимость Р от концентра- 
ции ядер ІР 35 в уран-ториевом сплаве для различ- 
ных значений начальной энергии нейтронов. 
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Рис. 3. Критерий возникновения неуправляемой цепной ре- 
акции в среде тп^пТН 232 в зависимости от ядерной 
концентрации урана, %: Р>1 - необходимое условие 
инициации; вероятность со рассчитана по формуле ( 8) 

Заключение 

Расчеты показывают, что во всех случаях в торие- 
вых сплавах допустимо присутствие ядер ІР 35 и Ри 239 в 
технологически значимых количествах. Даже в прене- 
брежении утечкой нейтронов (бесконечная среда) 
массовая концентрация ядер в сплаве может соста- 
влять 7 %, не создавая проблем ядерной безопасности. 

В компактных одиночных слитках в форме па- 
раллелепипеда с размерами 80x140x35 мм концен- 
трация делящихся ядер может быть увеличена до 
18 %, что делает их высококачественным полуфаб- 
рикатом для изготовления ядерного топлива. 

Полученные значения “времен жизни” нейтро- 
нов для различных видов размножающих сред 
играют важную роль в физике бустеров (подкрити- 
ческих сборок) и могут использоваться в процессе 
расчета а также эксплуатации указанных систем, в 
которых в качестве топлива применяются уран-то- 
риевые и плутоний-ториевые композиции. 
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